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1.0 SUHMARY 


This design note presents 2-lMU FOI performance data for the 
sequential probability ratio test (SPRT) durinp shuttle entry. Also 
included are current modelling constants and failure thresholds for 
the full mission 3B entry through landino trajectory. FDI results 
are presented in a “raw data" tabular format in order to furnish the 
reader with as much data tracking test detection and isolation per- 
formance data as is possible, with a minimal amount of data procession. 
Minimum 1002 detection/isolation failure levels and a discussion of 
the effects of failure direction are presented. Finally, a limited 
comparison of failures introduced at trajectory initiation sho\/s that 
the SPRT algorithm performs slightly worse than the data tracking test 
(Reference 1). 

2.0 IHTRODUCTIOfl 

Last September the SPRT algorithm was baselined at the Level F? 

OFT Entry SDR to perform the onboard 2 and 3 I?FU FDI testing with 
skewed IMU's. In order to both develop and verify the method, a 
subroutine incorporating the 2-IMll SPPT vias added to the PHIFDI triple 
string IMU simulation program on the JSC Univac 1110. This report 
contains an evaluation of the pres.*,..t SPP.T formulation (Reference 2) 
in detecting and identifying soft IflU failures. 

3.0 DISCUSSION 

The 2-IMU SPRT performance data presented in section ^.0 were 
generated by the IMUFDI program, version 18. C, interfaced with the 
2-IMU SPRT subroutine, as described in Reference 2. The follov/ing 
paragraphs contain error modelling data, guidelines, and other 
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constants specifying ^he exact conditions under which the enclosed 
failure test case results were generated. 

3.1 Guidelines 

• Refer^ice Mission 36 entry 

t Simulation begins at entry Interface (400 Kft.), 

and ends at touchdown* 1945 sec. after entry Interface 

I Only IMU #1 and IMU #2 are ON; IMU #3 Is downnoded throuohout 
all test cases 

t All failures are Introduced Into IMU #1 at tine t*0 (400 Kft.) 

t First detectlon/lsolatlon tests are perfomed at tine t»10 
sec., subsequent tests are perfomed every 5 sec. thereafter 

t Each failure case Is tested through 30. Monte Carlo cycles 

• The 11** nav base pitch Is modeled 

• The Kearfott IW glnbal sequence Is used (2YX-1nner middle 
outer) 

The following limitations of the IWJFOl program should be mentioned: 

• All IMU* s are assumed collocated. 

• No tangential or centripetal forces are modeled. 

• A 3 gimballed IMU error model Is used, error in the 4th 
(Inner roll) gimbal Is unnodeled. 

• This Is an open loop simulation 

3.2 Filter Constants 

The first order v/hitening filter is characterized by the follov»inn 


constants: 
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Autocorrelatlon time 
"'^GYRO * 

’Aca * 

Gains 

‘^YRO * 

K^Ccl * *0®^ 


3.3 Base Failure Thresholds 

The base failure thresholds are plotted In Figures 1 and 2» 
together with the 100 Monte Carlo cycle envelopes of nominal data 
before filtering. The gyro threshold Is a 3rd order polynomial 
function of time, specified by the following coefficients: 

TGYROq • 3.7 . E-4 

TGYRO^ » 3.4967 E -6 

TGYROg » -1.1786 E-10 

T 6 YRO 3 « -5.8263 


E-13 




TIac* see. frtn tntry teterfece 



lOOtOOO ft. W YeucMoi.li 

Yiee* see. fron Cntry tetcrfece 

i Figure 2. Incremental AV Differences 
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The accelerometer threshold Is a 3rd order polynomial which latches to 
a constant level at t«1145 sec. Before 1145, the threshold Is specified 
by the following coefficients: 


TACCLq » 

2.83 

E-5 

TACCL^ » 

4.3503 

E-8 

TACCLg * 

-5.3665 

E-11 

TACCLj • 

9.7743 

E-14 


3.4 log Likelihood Ratio Constants 

The failure thresholds on the residuals after filtering are also 
plotted In Figures 1 and 2, together with envelopes of 100 Monte Carlo 
cycle nominal residuals. These thresholds are calculated from the base 
failure thresholds as described In Reference 1 using the follovnng 
constants: 

Attitude Transient Percentage 

*"nb, gyro ' 

ACCL * 

Attitude Transient Detect Le* el 

^GIN * radians 
Residual Standard Deviation 

®GYR0 * 

Oaccl * ^*21-5 Km/sec 
Mean False Alarm Rate 

ALPHA » At/T = 5/5000 « 10‘^ 

Where T » 5000 sec., mean time between false alarms 
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3.5 Ske< Matrix 


The Ideal skewed transformation from IMU fl to UUJ #2 stable 


platform coordinate frames Is given by: 



-.5000000000 

-.8090169944 

.3090169944 


.80901 6994A 
-.3090169944 
.5000000000 


-.3090169044 

.5000000000 

.8090169944 


3.6 IMU Platfonn to Hav Base Euler Matrix 

The Kearfott stable platform to navigation base transformation 
matrix Is given by: 

C^Cd Ct|>S6 -S4» 

-c<i>se-t'S(t»SiiiCe S(t>C4» 

S(^S0+C(>S<{)C0 -S<^C0+C<^S^(>S0 

where S*s1ne> C*cos1ne» 0 are the X, Y, Z gimbal angles* 
respectively. 
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3.7 IMU Entry Error Model 



Engineering Values 

Program Values 

AC££LCnOK£TER ERRORS (la) per axis 



bias 

SO ug 

.490333369-006 Km 

scale factor 

100 pm 

.0001 

input axis misalignment 

IS arc sec 

.727220522-004 rad. 

quantization 

lx10“® km/sec 

1x10”^ Km/sec 

eVRO ERRORS (1o) per axis 



bias drift 

0.035 deg/hr 

.169684788-006 rad 



sec 

g-sensltive drift-input axis 

0.025 deg/hr/g 

.123592917-004 (rad/sec) 
(IQn/sec^) 

^-sensitive drift-spin axis 

0.02S deg/hr/g 

.123592917-004 (rad/sec) 
(Kta/sec^) 

2 

9 -sensitive drift-input/spin axis 

0.025 deg/hr/g^ 

.125826815-002 (rad/secl 
(Kw/sec^)^ 

scale factor 

200 PPM 

.0002 

counting alignment 

60 arc sec 

.290888209-003 rad 

IKU ERRORS (1a) | 



IKU to nav. base 

• 42.4 arc sec . 

.2056890249-003 rad 

gimbal non-orthoconalUy 

50 arc sec 

.2424068405-003 rad 

resolver bias term 

0. 

0 

resolver sinusoidal term 

30 arc sec 

.145^441043-003 rad 

resolver multiplicative speed 

2 

2 

gimbal quantization 

20 arc sec 

.969^73622-004 rad 

■ IMTIAL ALICraSNT ERRORS (lo) 


• 

each axis 

132 arc sec 

.6399540589-003 rad 


ORIGINAL PAGE IS 
OP POOR QUALirX 
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4.0 RESULTS 

Detection and Isolation sensitivities of the 2-IMU SPRT algorithm 
were tested on failures In gyro drift, accelerometer bias, and accelero- 
meter scale factor, Introduced Into IMU #1 at t®0 (400,000 ft). Failure 
levels were selected to correspond with the cases tested in Reference 2 
so that performance of the SPRT and data tracking test could be comparod 
under similar conditions. Failures In the following orientations were 
examined: 

Single axis: X, Y, and Z axes 
Dual axis: -45® and +45* In XY plane 
Finally, for the sake of completeness, SPRT performance In the presence 
of no failure Is summarized. 

Each failure case is summarized in a six column table. The first 
column Is the Monte Carlo cycle number. The second Is the time (sec. 

\ 

from entry Interface) of the first detection; this number will be equal 
to zero If there has been no detection during the cycle. The third 
column is the type of detection, ACCL or (lYRO; this field will be blank 
If there has been no detection. Columns 4 and 5 are the same as columns 
2 and 3, except that they pertain to the first isolation instead of the 
first detection. The sixth column is the IMU configuration control flag 
AFAIL, which is set after the first isolation. This number should be 
equal to 5 for all IMU ^1 failures; it would be equal to F if the isola- 
tion logic indicated an IMU 42 failure. 

4.1 Qyro Drift 

Tables 1, 2, and 3 contain FOI performance summaries for .5®/hr, 
l®/hr, and 2*/hr failures, respectively. The follovHng observations 


are made on the basis of these data: 
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• ,5*/hr detection vias high but never inof , ranging from 83?5 
In the Z axis case to 97?5 In the -45* case 

• l®/hr Isolation was IW, In only the case, falling to 
97% In the Y and Z cases, to 43% In the X case and 3% In the 
-45® case 

f 2®/hr Isolation was 100% In all but the -45® case. Average 
Isolation times were sharply reduced from the l®/hr cases. 

For example, Y axis average Isolation time fell from 1515 sec. 
to 479 sec. 

4.2 Accelerometer Bias 

Tables 4, 5, and 6 contain FOI performance summaries for lOOOug, 
IZOOpt-;, and 2000pg failures, respectively. The following observations 
are made on the basis of these data: 

t lOOOpg detection v»as high but never inr , ranglno from R7% 

In the -45® and Y cases to 93% In the X axis case 

• 1200pg detection vrais 100% In all cases. 1200pg isolation 
ranged from 0% In the -45® and X cases to 50% In the +45® case 

9 2000pg Isolation v/as 100% In the +45®, Y, and Z cases, falling 
to 20% In the X case and 3% In the -45® case. 

4.3 Accelerometer Scale Factor Error 

Tables 7, 8, and 9 contain FDI performance sunmarles for 3000ppm, 
SOOOppm, and SOOOppm failures, respectively. The following observations 
are made on the basis of these data: 

9 The 100% detection level lies betv/een 3000ppm and SOOOppm 

9 SOOOppm isolation was at best 00% In the X axis case 
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• SOOOppm Isolation was 1002 in the -45®, Y, and Z cases, fall inn 
to 97% in the +45* case, and to 632 in the X aj^is case. 

4.4 Honinal 

Table 10 contains the FOI performance surmary in the nominal case. 

In a separate run (not included) the algorithm was tested during 100 
nominal Monte Carlo cycles with no false detection or isolation. 

5.0 conclusions 

The following conclusions are made on the basis of the data contained 
in section 4.0. 

• For all cases run with the SPRT, there were no false detections 
in a 100 Monte Carlo cycle nominal case, and no incorrect 
isolations in the failure cases. 

• Failure detection v/as sensitive, for all failure orientations, 
with the follov/ing 1002 detection levels 

gyro drift : .6®/hr 

accelerometer bias - llOOug 

accelerometer scale factor ~ flOOOppm 

I For the IMU pair tested (#1, /*2), the best isolation geometry 
holds for the +45®, Y, and Z axis cases. Best geometry 1002 
isolation levels were found to be: 
gyro drift l®/hr 

accelerometer bias 1400pg 

accelerometer scale factor 6000ppm 

• -45® and X axis failures i i gyro drift and accelerometer Mas 
exhibited poor isolation K^rometry, since they were closest to 
the ambiguity line at -31.7® in the XY plane discussed in the 
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appendix of Reference 3. In the -45® case, l®/hr Isolation fell 
to 3X, and 20nopg Isolation fell to 3?!. 

Reference 1 presented the following lOOf, levels for the tracking test: 
Detection: .5®/hr, SOOyg, 40n0ppm 

Isolation: l*/hr, IZOOpg, SOOOppm 

Comparing these levels with those above for SPRT, the tracking test is, 
slightly more sensitive than the SPRT. A comparison of response tines 
shov/s the tracking test to be slightly faster In detectlon/lsolatlon 
than the SPRT. 
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Table 2. 1 ®/hr Gyro Drift, IMU #1 
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Table 3. 2 ®/hr Gyro Drift, IMU #1 
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Table 4. 1000 ug Accelerometer Bias, IMU #1 
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Table 5. 1200 vg Accelerometer Bias, IflU #1 
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Table 6. 2000 yg Accelerometer Bias, IMU #1 
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Table 7, 3000 ppm Scale Factor Error, IMU #1 
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Table 8. 5000 ppm Accelerometer Scale 
Factor Error, IMU #1 
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Table 9. 8000 ppm Accelerometer Scale 
Factor Error, IMU #1 
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Table 10. Nominal IflU’s #1 and #2 
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